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' Abstract. Most of the work done in the past on the integrability structure of the Classical 

D . Heisenberg Spin Chain (CHSC) has been devoted to studying the su{2) case, both at the 

continuous and at the discrete level. In this paper we address the problem of constructing 
. integrable generalized "Spin Chains" models, where the relevant field variable is represented 

' by a A'^xA^ matrix whose eigenvalues are the 7V"^ roots of unity. To the best of our knowledge, 

such an extension has never been systematically pursued. In this paper, at first we obtain 
the continuous N x N generalization of the CHSC through the reduction technique for 
Poisson-Nijenhuis manifolds, and exhibit some explicit, and hopefully interesting, examples 
' for 3x3 and 4x4 matrices; then, we discuss the much more difficult discrete case, where 

a few partial new results are derived and a conjecture is made for the general case. 
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■ 1 Introduction 
O . 

The Heisenberg ferromagnet has been one of the first integrable models investigated and solved in 
pH I the framework of the Inverse Scattering Method, both at a classical [H [21 [3] and at a quantum 

level. In particular, at the quantum level, through the Quantum Inverse Scattering Method 
[H [5], the original Bethe solution [6] has been rediscovered and generalized, and the solu- 
^ I tion of the isotropic one-dimensional spin chain was just the starting example in a list of 

■ a number of outstanding achievements obtained in this field in the 80's and in the 90's of 
the last century. The scope of the present contribution is twofold. On one hand, in Sec- 
tion 2 the geometric approach introduced in the field of infinite dimensional integrable sys- 
tems by Magri and coll. [TJ [U [8] will be used to extend to N x N matrices the old results 
on the SU{2) classical Heisenberg chain obtained about 20 years ago; accordingly, new hier- 
archies of integrable evolution equations will be derived, and a few explicit examples will be 
shown, which might have some relevance for studying classical analogs of the so-called Potts 
model \10\ \TT\ I12j . On the other hand, in Section 3 the intriguing issues concerning the 
lattice version of the model O [13] will be reexamined, a few additional results in the 2x2 
case will be given, and finally a conjecture concerning the general N x N case will be formu- 
lated. 



*This paper is a contribution to the Proceedings of the Workshop on Geometric Aspects of Integ- 
rable Systems (July 17-19, 2006, University of Coimbra, Portugal). The full collection is available at 
|http: / /www. emis.de/journals / SIGMA/Coimbra2006.html , 
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2 The classical Heisenberg chain on the line for N X N matrices 

Let us first recall the basic notions and results underlying the integrability of the classical 
Heisenberg ferromagnet in the usual su{2) setting. As a starting point, we will adopt the 
simplest point of view, where an integrable (matrix) system in 1 + 1 continuous dimensions is 
regarded as the compatibility condition for two linear differential (matrix) systems, defining the 
so called associated Lax pair: 

V'x = (2.1) 
^Pt = Vi;. (2.2) 

In (|2.ip . (j2.2|) ip, S, V are n x n matrix valued functions of the {x,t) coordinates, S is the 
field variable appearing linearly in the so-called principal spectral problem, V has an additional 
(rational) dependence on a complex parameter A, denoted as the spectral parameter. It would be 
natural to look at ip as an element of a Lie group, and at S, V as elements of the associated Lie 
algebra. More correctly, due to the extra dependence on A, we should speak of "loop" Groups 
and Algebras. However, in the sequel we will always assume S to belong to the open set of 
the invertible elements of the unital associative algebra of the N x N matrices with (generally) 
distinct eigenvalues. 

The compatibility condition of the above linear system ()2.ip . ()2.2p reads: 

St = jV. + [V,S] (2.3) 
which is also known as the zero-curvature condition for the connection given by the one form: 
a = XSdx — Vdt. 



Actually, in the very simple example under scrutiny, the compatibility condition (j2.3p amounts 
to the compatibility (i.e. vanishing of the Schouten bracket) of the two Poisson tensors: 



VV = [V,Sl QV = V, 



X 



which then define the bi-Hamiltonian pair "P, Q. 

As is well known from [3 O [H], the Poisson tensors V, Q are linear operators from the 
cotangent bundle T*M to the tangent bundle TM to our configuration manifold, which is just 
an affine hyperplane to the linear space of matrix-valued C°° functions of the real variable x 
fulfilling homogeneous boundary conditions. The duality form between TM and T*M is the 
usual trace form ((;/), a) := f^^Ty:{(j)a)dx. As Q is invertible, we can introduce on TM the 
tensor (Recursion operator): 



S, / 0(:= St) 



(2.4) 



that has zero Nijenhuis torsion (whence the name "Nijhenhuis tensor" |14j ) and gives rise to 
a hierarchy of commuting flows. Under the above conditions we will say that M is equipped 
with a VM structure. 

In order to construct a hierarchy of integrable evolution equations, we have to choose a sub- 
manifold M' C M such that the Nijenhuis tensor has maximal rank on it. As shown in [8] we 
can achieve this goal by going on restricting the action of TV to smaller and smaller submani- 
folds as far as the structure is preserved. In particular, it is always possible to restrict Af on 
a characteristic leaf of Im Q, and then again on a leaf of the distribution generated by ImAf^ , 
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where r (assumed to be finite) is the Riesz index of A/" [15]. As in our case Q is invertible, the 
only relevant restriction is on 

ImTV = lmV{(f> : (f) = [S,C(] for some a}. (2-5) 

It can be readily seen that ImM = luiM^, KevM = KeicAf^, VA; > 1. Hence r = 1 and then we 
can restrict the VM structure on a characteristic leaf of the distribution (12.51) : 



TM' = {(P : Tr = 0} (A: = 0, . . . , iV - 1) 
yielding the symplectic leaves: 

Tr5'= = 4. (2.6) 

The hierarchy obtained by acting with TV (j2.4p on the starting symmetry [5,(7], C being 
an arbitrary constant matrix, (the so-called "rotational flow") consists however of nonlocal 
equations. As it was already shown in [8] for 2x2 matrices, the way to get local equations is 
then to invert M, which is obviously possible thanks to the above restriction. 

In the following, for the sake of simplicity, we will simply denote by M^~^^ the inversion of 
the restriction of J\f on ()2.6p . 

The most important case, and the best studied one, is when S € su(2), namely it is a 2 x 2 
Hermitean traceless matrix: 

S = S-a, S-S = l, a = (cri,cJ2,cT3). 

Finding J\f~^ is a rather easy task. Indeed, from 



one gets 



S, 



^[S, + (cqI + ciS)x, 



where the second term in the above formula comes from Ker [S, ■]. Whence, imposing (j) to be 
tangent to the symplectic leaves (|2.6p . we derive: 



Co = 0, ci = --^ / tr [S,Sx 



so that: 

^ :=M(-^)^' = i [5, - ^5 I\t[S,S.]<P. (2.7) 

From (j2.7p one gets a local hierarchy [8], whose first members are: 
(t>o = 0, 0-1 = Sx, 4>-2 = [S, Sxx]-, — 
In particular the HSC Hamiltonian flow is 4>-2 = [S,Sxx\, namely: 

St = [S, Sxx] = [S,Vn\ ^ ^ = / T> {Sxf. (2.8) 

The novel result contained in this section is the explicit form of M^~^\ and of the associated 
local evolution equations and Hamiltonian functionals for N x N matrices s.t. = I, which 
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we have denoted as HSC at the N^^ roots of unity. The derivation is quite simple. Indeed, from 
the very definition of the Nijhenhuis tensor (|2.4|) one easily gets: 

\ 1=1 1=0 / 

The conditions: 
Tr5'=0 = O 

that ensure (j) to be tangent to (j2.6p . allow to determine through quadratures, up to constant 
integration factors: 

«=i 

For a given A^, the hierarchy of evolution equations splits into — 1 sub-hierarchies, with 
starting symmetries {S^)x {k = 1, . . . , N — 1). In the following, we provide the explicit form of 
the first members of the hierarchy for = 3 and N = 4. 

Example 1. = 3. There are two families of local flows, with starting points: 

4^ = S,, ct>f^ = {Sxf 
whence the first few members of simpler family read: 

The Hamiltonian TL^2^ reads (up to constant factors): 

J —oo 

Example 2. A^ = 4. There are three families of local flows, with starting points: 

whence the next flow in the simplest family reads: 

=3[S3,5i + [5,S5,S],, 
while the associated Hamiltonian functional TL^2^ is given by: 



J —oo 



Remark 1. The flow (|2.8p is obviously consistent with the restriction = I ^ N. However, it 
belongs to the integrable hierarchies only for N = 2. In fact, the Lie derivative of the Nijenhuis 
tensor ()2.4p in the direction of ^ := [S, Sxx] on an arbitrary vector (j) can be calculated from the 
formula: 

C4M(b) = {C^{M){<p) +MC^4> 
whence, with the proper replacements: 

(/:^(A/-)(,^)=3[5,[0,5,]]. (2.9) 
From (j2.9p we readily see that the Lie derivative of the Nijenhuis tensor in the direction (12. Sp 



vanishes only for the subhierarchies starting with (f)[^^ = Sx- Hence, for A^ larger than 2, the 
flow (12. sp . though being Hamiltonian and possessing infinitely many conserved quantities (and 
commuting symmetries) is certainly not completely integrable. 
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3 The discrete case: some results and a conjecture 

In contrast with the continuous case, we cannot give any new result for arbitrary N xN matrices. 
However, before focussing on the well known 2x2 case, we will write down the compatibility 
scheme and recall the abstract VAf structure in the differential-difference setting, as it has been 
elucidated in the past by various authors \13 \ \n \ \TE[. The two relevant linear problems read: 

1pm+l = {I + ><Sm)ll^m, 

'(Pm^t = VmiX)'4>m- (3.1) 

In (j3.ip tl^, S, V are N x N matrices, depending on a discrete variable m G Z, on a contin- 
uous variable t G M and on a complex parameter A. It has been shown [13^ \T7\ [18] that the 
compatibility of the above linear problems yields the VM structure: 

(j) := Va = [S, a], 

</.':=iV0=^[5(Rc/>)] + ^{5</.}, (3.2) 
where: 

£ + 1 

R- = £■ _ (^<^)m := 0m+l- 

We note that R is the so-called ILW operator, satisfying the modified Yang-Baxter equation in 
the algebra of shift operators [T7| . 

Remark 2. Notice that now ImTV is larger than ImV; however Im'P is still an invariant 
submanifold for ImA/", namely Tr S^4> = implies Tr S'-Mcj) = V/. Indeed, from (|3.2|) it follows 
the recurrence relation Tr SWc/) = Tr 5'+ V, entailing Tr [M^ - 2:](/> = = Tr 5' [Af^ -I](t),yi. 

As we mentioned before, to the best of our knowledge explicit evolution equations can be 
found in the literature only for 2x2 matrices. They are for instance reported in the fundamental 
reference [2], where the corresponding "Lax pair" ()3.1I) is also shown. As it was first pointed out 
in [13], both the Discrete Heisenberg (DH) ferromagnet and the Discrete Translational (DT) flow 
belong to the kernel oiM^ — T, which is in fact the linear span of them. They read respectively: 

Sm,T — Pm{Sm + 'S'm+l) ~ Prn—l{Sm + 'S'm— l)) 

where Pm '■= tv{i+s s +i) • ^^^^ the above two equations can be derived from the other by 
applying the recursion operator M. 

The above results, and Remark 2 strongly suggest the following 

Conjecture. For N x N matrices obeying the cyclic constraint = I, the local evolution 
equations span Ker (A/"^ ~^)/ hence, starting from any of them, the others can be obtained by 
applying subsequently M. In other words, we can always find in Ker {J\f^ —2) a cyclic basis (f)^^^ 
such that: 0^'^+^) = Af^^''\ (f>^''+^^ = (t>^'^\ 
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